Using GENERATINGFUNCTIONOLOGY to Enumerate Distinct-Multiplicity Partitions 
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In fond memory of Guru Herbert Saul WILF (28 Sivan 5691- 12 Tevet 5772) zecher gaon Vbracha 



About a year ago, Herb WilffWl] posed, on-line, eight intriguing problems. I don't know the answer 
to any of them, but I will say something about the sixth question. 

Herb Wilf 6th Question: Let T(n) be the set of partitions of n for which the (nonzero) multi- 
plicities of its parts are all different, and write f(n) = \T(n)\. See Sloane's sequence A098859 for 
a table of values. Find any interesting theorems about f(n) . . . 

First, I will explain how to compute the first few terms of f(n). Shalosh can easily get the first 
250 terms, but as n gets larger it gets harder and harder to compute, unlike its unrestricted cousin, 
pin). I conjecture that the fastest algorithm takes exponential time, but I have no idea how to 
prove that claim. I am impressed that, according to Sloane, Maciej Ireneusz Wilczynsk computed 
508 terms. 

Recall that the generating function for the number of integer partitions of n whose largest part is 
< m, p m (n), is the very simple rational function 



The main purpose of this note is to describe, using Generatingfunctionology, so vividly and lucidly 
preached in W's classic book [W2], how to compute the generating function (that also turns out 
to be rational) for the number of partitions of n whose largest part is < m and all its (nonzero) 
multiplicities are distinct, let's call it f m (n). As m gets larger, the formulas get more and more 
complicated, but we sure do have an answer, in the sense of the classic article [W3], for any fixed 
m, but of course not for a symbolic m. 

Even more is true! Because, like (1 _ g)(1 _ g 2).,. (1 _ g ^) , the generating function of / m (n), Y^=o fm(n)q n , 
turns out (as we will see) to only have roots-of-unity poles, whose highest order is m, it follows 
that fm(n) is a quasi-polynomial of degree m — 1 in n. Now that's a very good answer! (in W's 
sense, albeit only for a fixed m). 
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Preamble 




How to Compute Many terms of /(n)? 

p m (n) is very easy to compute. For example, one may use the recurrence 

[n/m\ 

Vm(n) = p m -i(n) + ^ p m -i(n-mi) , 

i=l 

together with the initial condition p\(n) = 1, p m (0) = 1. 

How can we adapt this in order to compute / m (n)? The contribution from the partitions counted 
by fmin) where m does not show up is / m _i(n), in analogy with the p m -i(n) term in the above 
recurrence. But if m does show up, it does so with a certain multiplicity, i, say, where 1 < i < 
[n/m\, and removing these i copies of m results in a partition counted by / m _i(n — mi) -so all 
its multiplicities are different- and in addition none of these multiplicities may be i. Continuing, 
we are forced to introduce a much more general discrete function f m (n; S) whose arguments are m 
and n and a set of "forbidden multiplicities" , S. 

So let's define f m ( n ', S) to be the number of partitions of n with parts < m, with all its multiplicities 
distinct and none of these multiplicities belonging to S. Our intermediate object of desire, f m (n), 
is simply / m (n;0), and the ultimate object, f(n), is / n (n;0). 

The recurrence for f m {n]S) is, naturally 

[_n/m\ 

f m (n;S) = f m - 1 (n;S)+ ^ f m -i{n - im; S U {i}) , 

i=l,i£S 

because once we decided on the number of times m shows up, let's call it i, where i is between 1 
and [m/m\ and i S, the partition (of n — mi) obtained by removing these i copies of m must 
forbid the set of multiplicities S (J {i}. 

In the Maple package DMP, procedure qnmS(n,m,S) implements f m (n;S) and procedure qn(n) 
implements f{n). 

Inclusion-Exclusion 

Let V m (n) be the set of partitions of n whose parts are all < m, in other words, the set that p m (n) 
is counting. Consider the set of all partitions whose largest part is < m, where we write a partition 
in frequency notation: 

V m := {l ai 2 a2 ...m a '" |ai,...,o m > 0} . 

For example 1 3 2 5 4 2 is the partition of twenty-one usually written as 4422222111. Introducing 
symbols x\, X2, ■ ■ ■ , x m , we define the Weight of a partition to be 

Weight(l ai 2 a2 ...m am ) := x^x^ 2 • • • <7 . 
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The weight- enumerator of V m is, by ordinary- generatingfunctionology 

Weight(V m ] ' 



(1 - - x 2 ) ■ ■ ■ (1 - x m ) 

since we make m independent decisions: 

how many copies of l?(Weight enumerator = 1 + x\ + x\ + . . . = (1 — , 

how many copies of 2? (Weight enumerator = l+x 2 + X2 + ... = (l — x 2 ) _1 ) , 



• how many copies of ml (Weight enumerator = 1 + x m + x 2 m + . . . = (1 — x m ) 1 ). 
But we want to find the weight-enumerator of the much-harder-to- weight-count set 

F m :={l ai 2 a2 ...m a ™\a 1 ,...,a m >0;a i ^a J {if a i >0,a j >0)} . 

Calling the members of T m good, we see that a member of V m is good if it does not belong to any 
of the following (™) sets, SV,- 1 < i < j < m: 

S i:j := {l ai 2 a2 . . . m am G V m \<n = a 3 > 0} . 

By inclusion-exclusion, the weight-enumerator of T m is 

Ei-^^ght^^S^ , 

G 

where the summation ranges over all 2 m ( m " 1 )/ 2 subsets of {{i,j)\l<i<j<m}. 



But the G's can be naturally viewed as labeled graphs on m vertices. Such a graph has sev- 
eral connected components, and together they naturally induce a set partition {C\, C 2 , . . . , C r } of 
{1,2, ...,m}. We have: 



Weight (^.. ^ ^ Si^j = Y[weight(Ci) , 



where if \S\ = 1, S = {s}, say, then weight(S) = j^r, and if |5| = d > 1, S = {s\, s 2 , . . . Sd}, say, 
then 

rv* ry* ■ ■ ■ / T i 

weight(S) = _ S1 32 . 

i x Sl x S2 • • • x Sd 

To justify the latter, note that if vertices s\, s 2 , ■ ■ ■ Sd all belong to the same connected component 
of our graph then, by transitivity, we have that all a Sl = a S2 = ... = ... a Sd > 0, and the weight- 
enumerator is the infinite geometric series 

oo 

E, \a X Sl X S2 ■ ■ ■ X Sd 

\ X Sl ' ' ' x s d ) 



a = l 



1 Xc-. X 
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But quite a few graphs correspond to any one set-partition. To find out the coefficients in front, 
for any set-partition {Ci, C2, • • • , C r } of {1, ... , m} we must find 

£(-i) |G1 , 

G 

summed over all the graphs that gives rise to the above set partition. But this is the product of 
the analogous sums where one focuses on one connected component at a time, and then multiplies 
everything together. 

Let's digress and figure out J2g(~ 1)' G over an connected labeled graphs on n vertices. For the sake 
of clarity, let's, more generally, figure out ^2 G y^ G ^ with a general variable y. 

By exponential-generatingfunctionology[W2] (see also [Z]), this sum is nothing but the coefficient 
of t n /n\ in 



log (jT(i + y)G): 



,i=0 



Going back to y = — 1, we see that we need the coefficient of t n /n\ in 

00 00 

= E(- 1 )"- 1 ^ = E(- i r-v-D^ . 

n=l n=l 

So the desired sum is (— l) n ~ 1 (n — 1)!. 

Let's define for any set of positive integers, S, 

mishkal(S) ~ { 1/(1 " Xs) ' if 151 = 1 Where S = {S} 5 

1 j " I (-1)^ - l)!(x Sl • • • x Sd )/(l - x Sl ■ ■ ■ x Si ), if |S| = d > 1 where S = {s u ... s d }. 

For any set partition C = {C\, . . . , C r } let's define 

Mishkal(C) = mishkal(Ci) ■ ■ ■ mishkal(C r ) 
It follows that the weight-enumerator of T m according to W eight{l ai 2° 2 . . . m am ) := x^x^ 2 • • • x 



rn 



IS 



^Mishkal(C) , 



where the sum has B m terms (B m being the Bell numbers), one for each set-partition of {1, ... , m}. 

Finally, to get an "explicit" formula (as a sum of B m terms, each a simple rational function of q), 
for the generating function Yl^Lo fm( n )l n , all we need is replace Xi by q l , for i = 1 . . . m, getting 

n=0 C 
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where for a set partition C = {Ci, . . . , C r } 



Poids(C) = poids(C\) ■ ■ • poids(C r ) 



and where for an individual set S: 



poids(S) = 



(-if-^d- l)!g Sl+ - +Sd /(l - q 



,Si + ... + S d 



if \S\ = 1 where S = {s} ; 
), if \S\ = d > 1 where S = {si, . . . s d }. 



It follows that indeed / m (n) is a quasi-polynomial of degree m — 1 in n. Furthermore, since the 
only pole that has multiplicity m is q = 1, it follows that the leading term (of degree m — 1) is a 
pure polynomial. 

The generating function, Y^=o fm{n)q n ■, for any desired positive integer m, is implemented in 
procedure GFmq(m,q) in the Maple package DMP. For the Weight-enumerator (or rather with Xj 
replaced by q l Xi, for i = 1, . . . , m), see GFmxq(m,x,q) . Since the Bell numbers grow very fast, the 
formulas get complicated rather fast, but in principle we do have a very nice answer for any specific 
m, but in practice, for large m it is only "nice" in principle. Of course it is anything but nice when 
viewed also as function of m, and that's why f(n) = f n (n) is probably very hard to compute for 
larger n. 

To see the outputs of GFmq(m,q) for 1 < m < 8 see: 

http : / / www . math . rutgers . edu/~zeilberg/tokhniot/ oDMP3 . 

Asymptotics 

Recall that Hardy and Ramanujan tell us that as n goes to infinity, pin) is asymptotic to ^^ exp(C^fn) 
where C = \j2/2>n = 2.565099661 . . ., and hence \ogp(ri) / y/n converges to C. By looking at the 
sequence log /(n)/y / n for 1 < n < 508, it seems that this too converges to a limit, that appears to 
be a bit larger than 1.517 (but of course way less than 2.565099661 . . .). Let's call that constant 
the Wilf constant. 

The numerical evidence is here: http://www.math.rutgers.edu/~zeilberg/tokhniot/oDMP4 . 
Let me conclude with two challenges. 

• Prove that the Wilf constant exists. 

• Determine the exact value of the Wilf constant (if it exists) in terms of ir or other famous 
constants. Failing this, find non-trivial rigorous lower and upper bounds. 
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